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Abstract 

Within the functional calculi of Bochner-Phillips and Hirsch, we describe the operators 
of distributed order differentiation and integration as functions of the classical differenti- 
ation and integration operators respectively. 

1. Introduction and Preliminaries. In the distributed order calculus [11], used in 
physics for modeling ultraslow diffusion and relaxation phenomena, we consider derivatives 
and integrals of distributed order. The definitions follows. 

Let /i be a continuous non-negative function on [0, 1]. The distributed order derivative D^^^ 
of weight for a function (p on [0, T] is 



where D^"^ is the Caputo-Dzhrbashyan regularized fractional derivative of order a, that is 



1 




(1) 







1 d 
r(l - a) It 



(t-r)-"(^(r) rfT-r"(^(0) 



< t < T. 



(2) 







Denote 



1 




(3) 
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It is obvious that A; is a positive decreasing function. The definition (l)-(2) can be rewritten as 

t 

(D(/^)y,) {t) = jj k{t - tMt) dr - k{tMO)- (4) 



The right-hand side of (4) makes sense for a continuous function ip, for which the derivative 

d * 

— / k{t — t)lp{t) dr exists. 
dt 

If a function (/? is absolutely continuous, then 

t 



Below we always assume that /i G C^[0, 1], //(I) ^ 0, and either /i(0) ^ 0, or ~ aa^, 
a, > 0, as a — > 0. Under these assumptions (see (llj). 

~ s~^(logs)^V(l), s^O, 

~ -s-2(logs)-V(l), s^O, 

so that G Li(0,T) and A; does not belong to any Lp, p > 1. We cannot differentiate under 
the integral in (4), since k' has a non-integrable singularity. 

It is instructive to give also the asymptotics of the Laplace transform 



lC{z) = I k{s)e-''ds. 

Using (4) we find that 



1 

a-l , 



}C{z) = J z"~ fi{a) da, 



so that IC{z) can be extended analytically to an analytic function on C \ ]R_ = {z E C : 
Imz = 0,Rez < 0}. If z G C\R_, \z\ ^ oo, then [nj 

/C(z) = f^ + 0((log|z|)-2); (6) 
log z 

see [11] for further properties of /C. 

The distributed order integral I^^^ is defined as the convolution operator 

t 

(t) = j x{t- s)f{s) ds, 0<t<T, (7) 
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where x(t) is the inverse Laplace transform of the function z 



zJCiz) 

7+ioo 

f '. 

dt2ni J z zlC{z) 



M^) = :^7r- / —■-l^.dz, 7>0. (8) 



7— lOO 

It was proved in [11] that x G C°°(0, 00), and x is completely monotone; for small values of t, 

x(t)<Clogi |x'(t)| <Ct-Mog^, (9) 

If / G Li(0, T), then f = f. 

The aim of this paper is to clarify the operator-theoretic meaning of the above constructions. 
It is well known that fractional derivatives and integrals can be interpreted as fractional powers 
of the differentiation and integration operators in various Banach spaces; see, for example, 

[D sua [15]. 

du 

Let A be the differential operator Au = in Lp{0,T), 1 < p < 00, with the boundary 

condition ■u(O) = 0. Its domain D{A) consists of absolutely continuous functions u G Lp{0,T) 
, such that 'u(O) = and u' G Lp(0,T). We show that on D{A) the distributed order differen- 
tiation coincides with the function C{—A) of the operator —A, where C{z) = z]C{z), and the 
function of an operator is understood in the sense of the Bochner-Phillips functional calculus 
(see [iniElIIS]). 

Moreover, if p = 2, then the distributed order integration operator I^^^ equals A/'(J), where 

1 * 
Aflx) = , J is the integration operator, {Ju){t) = J u{T)dT. This result is obtained 
C{x) 

within Hirsch's functional calculus [71 [8] giving more detailed results for a more narrow class 

of functions. As by-products, we obtain an estimate of the semigroup generated by —C{—A), 

and an expression for the resolvent of the operator I'^'^^. 

2. Functions of the differentiation operator. The semigroup Ut of operators on the 
Banach space X = Lp(0,T) generated by the operator A has the form 



(Utf) (x) 



f{x-t), iiO <t<x<T; 
0, if < X < t. 



X G (0,T), f > 0. This follows from the easily verified formula for the resolvent R{X,A) = 
{A — A/)~^ of the operator A: 

X 

{R{X, A)u){x) = - [ e-'^^-y^u{y) dy; (10) 



see [10] for a similar reasoning for operators on Lp(0,oo). The semigroup Ut is nilpotent, 
Ut = for t > T; compare Sect. 19.4 in [6j. It follows from the expression (10) and the Young 
inequality that ||-R(A,v4)|| < A > 0, so that Ut is a Co-semigroup of contractions. 
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In the Bochner-Phillips functional calculus, for the operator A, as a generator of a contrac- 
tion semigroup, and any function / of the form 



f{x) = J {1- e~*'')a{dt) +a + bx, a,b>0, (11) 



where cr is a measure on (0, oo), such that 

oo 


the subordinate Co-semigroup U/ is defined by the Bochner integral 

oo 

Ui = J {Usu)at{ds) 



where the measures at are defined by their Laplace transforms, 

oo 



The class B of functions (11) coincides with the class of Bernstein functions, that is functions 

/ G C([0, oo) n C°°(0, oo), for which /' is completely monotone. Below we show that £ G i3. 
The generator of the semigroup [// is identified with —f[—A). On the domain D[A), 

oo 

A-^u^ -au + bAu + J{UtU-u)a{dt), u e D{A). (12) 



Theorem 1. (i) Ifue D{A), then A^u = -D(^)m. 

(ii) The semigroup Uf decays at infinity faster than any exponential function: 

\\U^\\< ae"^* for any r > 0. (13) 

The operator A'' has no spectrum. 

(iii) The resolvent R{X, —A^) of the operator —A^ has the form 

X 

iRiK-Ac^i.)^jn(.-.Hs)as. uex, (14) 



where ^ ^ 

rx{s) = );^«a(s), (15) 



and u\ is the solution of the Cauchy problem 

B>^^\^ = Xu^, u^{0) = 1. (16) 

(iv) The inverse [—A'-') coincides with the distributed order integration operator 

(v) The resolvent of I*^^^ has the form 

(l(^) - Xiy^u = -^u - -iri/, A ^ 0. (17) 

Proof Let a{dt) = -k'{t)dt. By (3), 

1 

f at~°'^^ 
k'it) = - ^ _ x ^(q) c?a, 



so that 

oo 1 oo 

t f afi(a) , f t~" , 

a(at) = / — r eta / at. 

1 + t ^ ^ y r(l-a) J l + t 



Using the integral formula 2.2.5.25 from [T3] we find that 



< oo. 



ft , , , f a^ia) , 

/ aidt) = TT I — — -da 

J l + t ^ ' J (sina7r)r(l -a) 



Let us compute the function (11) with a = b = 0. We have 

oo oo 

fix) = - j (1 - e^*^) k\t) dt = x j e-'^'kit) dt = xlC{x) = C{x). 



The corresponding expression (12) for A'-'u, u G D{A), is as follows: 

oo X 

[A'^u) {x) = - I [{Utu){x) - u{x)]k\t) dt = - [u{x - t) - u{x)]k'{t) dt + u{x) I k'{t) dt 



= —k{x)u{x) — J [u{x — t) — u{x)]k'{t) dt. 



By (4), we find that A^u = -ID^'^'^u, u G D(A). 

The function Ci^z) is holomorphic for Re 2; > 0. We will need a detailed information (refining 
(6)) on the behavior of Re£(a" + «r), a, r G M, a > 0, when |r| 00. We have 



Re£(cr + zr) = ^ ip{a,a,T)^{a) da 
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where 



^{a, a, t) = ((T^ + r^)"/^ cos (^a arctan — j . 
We check directly that ^p{a, a, 0) = cr", 

dip{a,a,T) 2 , 2^a/2-l ( ^ ^ ( + ^ 

— = aia + T ) ' cos a arctan — r — a tan a arctan — 

OT \ a/ I \ a 



> 0, 



d { (X (T t1 

and — — > for a < 1. This means that the function Qair) = Re£(o" + it) (which is 

even in r) is strictly monotone increasing in r for r > 0. Its minimal value is 

1 

5'a(0) = j a'^nia) da. 



On the other hand, 

1 n/2 

ReC{a + iT)> y (a^ + r')"/^ cos da = ^ j {a^ + ry^'' fi (^^^ cost (it 



7r/2 tt/2 

= - [ e'^'fi ( costdt = -e^"/2 [ e-^'fi (l~-s] sin sds 

n J \7r J TT J \ TT J 



where q = ^ log(c"^ + t^)- By Watson's asymptotic lemma (see [I2]), since /i(l) ^ 0, we have 

tt/2 

e^'^^fi I 1 s I sins (is ~ Cq~^ 

V ^ J 

where C does not depend on a, r. Roughening the estimate a little we find that 



, 1 



e 



where < £ < | can be taken arbitrarily, and the positive constants C and p do not depend 
on a and r. 

It follows from (18) (see f3]) that for each t > the function x ^— e~^^^^^ is represented by 
an absolutely convergent Laplace integral. This means that the measure at{ds) has a density 
m(t, s) with respect to the Lebesgue measure. Moreover, 



7-t-/o<j 

i(t, s) = — / e"^ ■ e-*^(") rfz, 7 > 0. (19) 
2vrz J 



7+ioo 

1 

m( 

' ' 27r' 

7— lOO 

Since ?7t = for t > T, we have 

T 

[//w = y (f/,^x)m(t, s) ds. (20) 
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The representation (19) yields the expression 



We have 



oo 


1 

£(7 + zr)= {-J + iT)"fi{a) da, < r < cxd. 



oo 

\m{t,s)\ <— j e-'s^'^^Ur. 



The above monotonicity property of makes it possible to apply to the last integral the 
Laplace asymptotic method [L2l. We obtain that, for large values of t, 

\m{t,s) \ < Cr^e^'e-*3^(°). 

Changing 7 and C we can make the coefficient 5^7(0) arbitrarily big. By (20), this leads to the 
estimate (13). 

Due to (13), the resolvent 

00 

R{X,A^) = - [ e-^'Ufdt, (21) 



is an entire function, so that has no spectrum. 
It follows from (21) that 



00 

R{X,-A^) = [ e^%^dt, 

and if G X, Re A < 0, then 

[R{X, —A^)uj (x) = / e'^*(it / u{x ~ s)m{t, s) ds = / rx{x — s)u{s) ds 



where 

rx{s) = I e^'m{t,s)dt. (22) 



00 



For a fixed u G (|, 1), let us deform the contour of integration in (19) from the vertical line 
to the contour S^^^ consisting of the arc 

T^,u; = {z E C : |2;| = 7, I argz| < ujtt}, 

and two rays 

^^,Lu = {z eC: I argz| = icjvr, \z\ > 7}. 
7 



The contour S^^^^ is oriented in the direction of growth of arg z. By Jordan's lemma, 



m{t, s) = — I e''- e-*^(") dz. 
2Tii 



Under this integral, we may integrate in t, as required in (22). We find that 



rx{s) 



1 



2Tn J C{z)-\ 

s. 



dz, s > 



(23) 



(for Re A > 0, 7 should be taken big enough). 

If A = 0, the right-hand side of (23) coincides with that of (8) (see also the formula (3.4) in 
[11]). and we prove that (-y4^) ^ = l(^). 

For A 7^ 0, we rewrite (23) as 



rx{s) 



2Tvi\ 



Ciz] 



C{z)-\ 



dz 



27rzA 



e^"" dz. 



(24) 



For < s < T, we have 



e^* dz = — lim 

i?,^oo 



e""" dz, 



\z\=R 
ti;7r<l arg2;l<7r 



dz 



\z\=R 
j-K< \ arg z|<7r 



7r 

< 2i?y" e^"'''''^dip < 2Rn{l 



Rs cos L 



as i? oo. 

Thus, the second integral in (24) equals zero, and it remains to compare (24) with the 
formula (2.15) of |TT] giving an integral representation of the function ux. 

The formula (17) follows from (15) and the general connection between the resolvents of an 
operator and its inverse ([lOj, Chapter 3, formula (6.18)). The theorem is proved. 

Note that the expression (17) for the resolvent of a distributed order integration operator is 
quite similar to the Hille-Tamarkin formula for the resolvent of a fractional integration operator 
(see |6], Sect. 23.16). In our case, the function ux is a counterpart of the function z i— > £'^(Az") 
(for the order a case). However, in our situation no analog of the entire function (the 
Mittag-Leffler function) has been identified so far. Accordingly, our proof of (17) is different 
from the reasoning in [6]. 

3. Functions of the integration operator. In this section we assume that p = 2. 

Hirsch's functional calculus deals with the class TZ of functions which are continuous on C \ 
(— oo, 0), holomorphic on C\ (— oo, 0], transform the upper half-plane into itself, and transform 
the semi-axis (0, oo) into itself. The class TZ is a subclass of B. 
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Another important class of functions is the class S of Stieltjes functions 

oo 

f(z) = a + I ^gC\(-oo,0], 



(ip(t) 

where a > 0, p is a non-decreasing right-continuous function, such that f < oo. If f is a 

- , H & & , i l+t 

nonzero function from S, then the function 



also belongs to S. 

If / G 5, then the function Hf{z) = f{z^^) belongs to TZ. It has the form 



oo 



Hf{z) = a+ I ^^rfp(A), 2:gC\(-oo,0]. 

For some classes of linear operators V, the function Hf{V) is defined as a closure of the 
operator 

oo 

Wx = ax + J V{I + XV)-^xdp{X), x E D{V). 



In particular, this definition makes sense if —V is a generator of a contraction Co-semigroup, 
and in this case the above construction is equivalent to the Bochner-Phillips functional calculus 
[21 H]. In addition, by Theorem 2 of [7j, if {—V)~^ is also a generator of a contraction Co- 
semigroup, then 

[Hj{V)]-' = H^{V-'). (25) 
In order to apply the above theory to our situation, note that [7] 

oo 

z'' = ——^ [^^X--dX, 0<a<l, 

r(a)r(l — a) J 1 + Xz 





whence 



where 



oo 




1 



Thus C{z) = Hf{z), with 



oo 
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It follows from Watson's lemma [I2j that P{X) < C(logA) ^ for large values of A. Therefore 

oo 

/ a\ < oo. 

J 1 + A 



Denote Af{z) = Hj{z) = 

\i V = —A, then {—V)^^ = —J, where J is the integration operator. It is easy to check 
that ((J + J*)u, > ((■,■) is the inner product in L2(0, T)). Therefore — J is a generator of 
a contraction semigroup. 

After these preparations, the equality (25) implies the following result. 

Theorem 2. The operator I^''^ of distributed order integration and the integration operator J 
are connected by the relation 
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